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Recent observations of oscillatory features in the optical response of photosynthetic complexes 
have revealed evidence for surprisingly long-lasting electronic coherences which can coexist with 
energy transport. These observations have ignited multidisciplinary interest in the possible role of 
quantum effects in biological systems, including the fundamental - though still unresolved - ques- 
tion of how electronic coherence can survive in biological surroundings. Here we show that in 
photosynthetic complexes, non-trivial spectral structures in environmental fluctuations can allow 
for non-equilibrium processes that lead to the spontaneous generation and sustenance of electronic 
coherence even at physiological temperatures. Developing advanced new simulation tools to treat 
these effects, this new insight provides a firm microscopic basis to successfully reproduce the experi- 
mentally observed coherence times in the Fenna-Matthews- Olson complex, and thus sets the ground 
for the future assessment of the role of quantum effects in photosynthetic light harvesting efficiency. 



Introduction - Photosynthesis is a fundamental biolog- 
ical process that provides the primary source of energy 
for almost all terrestrial lifffl In its early stages, ambi- 
ent photons are absorbed by optically active molecules 
(pigments) in an antenna complex, leading to the forma- 
tion of molecular excited states (excitons). These then 
migrate by excitation energy transfer (EET) through 
pigment-protein complexes (PPCs) to a reaction centre 
(RC) where the exciton's energy is used to release an elec- 
tron - see Fig. [l] Remarkably, these processes often have 
a quantum efficiency of almost 100%^ ^, and uncovering 
the underlying biological design principles could inspire 
important new developments in artificial light harvest- 
ing technologies^. The potential novelty of a biomimetic 
approach to light-harvesting is underlined by the recent 
and unexpected observation of robust, long-lasting oscil- 
latory features in 2-D spectra of PPCs extracted from 
bacteria, algae and higher plants. Using ultrafast non- 
linear spectroscopy, sustained beating between optically 
excited states lasting several hundreds of fs at room tem- 
perature, and up to nearly 2ps in the Fenna-Matthews- 
Olson (FMO) complex at 77K, have been observecP^^. 
These experiments have been interpreted as evidence 
for electronic coherences between excitons, with lifetimes 
which are, surprisingly, over an order of magnitude longer 
than coherences between electronic ground and excited 
states^. Such coherence times are long enough for EET 
and excitonic coherence to coexist in these PPCs, condi- 
tions under which a sophisticated interplay of quantum 
and dissipative processes theoretically optimises trans- 
port efficienc}EHi^. 

While many proposals for how quantum effects might 
enhance biological light-harvesting have been advanced 
over the last five years, most of these have used simple, 
phenomenological methods to include decoherence^^^. 
In this letter we tackle the more fundamental prob- 
lem of elucidating the microscopic mechanisms through 
which the wide range of electronic coherence times ob- 
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FIG. 1: The generic organisation of the early stages of light- 
harvesting in natural photosynthesis. Excitons created in the 
antenna complexes migrate via dipolar coupling between chro- 
mophores in different pigment-protein complexes and are fi- 
nally transferred to a reaction centre where charge separation 
occurs. Under the low light conditions where light-harvesting 
is most efficient, most of the intermediate pigment-protein 
complexes transport only single excitons at a time and each 
complex functions independently. Chromophores are modeled 
as having a ground state \g) and optically excited state \i). 
Couplings between chromophores leads to the formation of de- 
localised excitonic eigenstates |ei), with transitions between 
these states mediated by environmental fluctuations 



served in PPCs actually arise. By including full tem- 
perature dependence into the numerically exact Time 
Evolving Density with Orthogonal Polynomial Algorithm 
(TEDOPA) technique of Refs.^"^ ^^, we show that the 
presence of resonant structure in the spectral function 
of PPC chromophores leads to non-equilibrium quantum 



dynamics in which picosecond electronic coherence can 
be driven and supported by quasicoherent interactions 
between excitons and spectrahy sharp vibrational envi- 
ronment modes^^H^. At the same time the model also 
predicts, for the same parameter values, that both the 
decay of ground-excited state coherences and the excita- 
tion energy transport exhibit the correct time scales, thus 
explaining three crucial observations within the same 
model. We complement our numerics with an approx- 
imate semi-classical model that illustrates the essential 
physical mechanisms by which discrete modes, driven 
far from equilibrium by exciton injection, may sponta- 
neously generate and sustain oscillatory EET and elec- 
tronic coherence against aggressive background decoher- 
ence, and which also helps to isolate the roles of other 
processes, such as non-Mar kovian pure dephasing^^ ^^. 
Under these conditions, energy transport proceeds in 
a non-equilibrium fashion, transiently violating the de- 
tailed balance conditions which are often invoked to set 
fundamental limits on light harvesting efficiency^. The 
conditions that allow for the observation of this striking 
behaviour in ensemble experiments are discussed. 

Model - Let us consider the following standard Hamil- 
tonian of a PPC consisting of a network of chromophores 
which are linearly coupled to environmental fluctuations 
in their excited states^ ^^. These fluctuations can be due 
to vibrations in the surrounding protein matrix and/or 
intrachromophore vibrational modes. Deflning bosonic 
annihilation and creation operators a^/c, aj^ for the k — th 
independent oscillator coupled to site i, and denoting an 
optical excitation on site i by the state |i), the Hamil- 
tonian {h = 1) for a single electronic excitation can be 
written as 

H = J2{{e.+X,)\i){i\+H,^)+J2j,,\i){j\, (1) 

i i^j 

where coupling to the vibrational modes of the environ- 
ment induces local energy fluctuations via the mode dis- 
placement operators Xi = Yjk9ik{ciik + «L) ^^^ ^F = 
^k^kCi\]^Ciik^ with ujk being the frequencies of the vi- 
brational fluctuations and Qik their coupling strength to 
state \i). The optical excited states \i) have coherent 
coupling matrix elements (arising from dipolar coupling) 
Jij = Jji. This coupling leads to delocalised excitonic 
eigenstates \en) with energies En (we order the states 
such that El is the lowest energy excitation) - see Fig. 
[l] Finally, we assume, for simplicity only, identical, in- 
dependent vibrational environments on each site charac- 
terised by a spectral density J{uj) = ^]^gij^S{uj — uok)- 
The key flndings presented here arise from our consider- 
ation of structured spectral functions which contain two 
contributions; a smooth background describing fluctua- 
tions, likely due to the protein environment, and also 
couplings to discrete vibrational modes, which may be of 
intramolecular origin. 

Physical mechanism of vibration-induced electronic co- 
/lerence-Rewriting the Hamiltonian in the excitonic basis 
of eigenstates {|en)}, so that |^) = Xl^ C^|en), it is shown 



in the Appendix | A 3| that the exciton-environment inter- 
action develops a transverse (non-adiabatic) term which 
couples different excitonic states through the bath dis- 
placement operators. Neglecting non-mar kovian effects, 
the background environment leads to incoherent exciton 
relaxation. However, coupling to quasi-resonant discrete 
modes allows for the possibility of coherent and reversible 
inter-exciton transitions via the long-lasting mechanical 
(coherent) motion of the mode displacements. 

The coupling to the discrete modes is described by the 
Hamiltonian terms. 



Hi = ^ ^{Qnm\en){em\ + h.c), 
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Qnm = / , v Sk^kCnCl^{aik + g^fc), 
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where uok are the frequencies of the coherent modes and 
Sk = (gk/^k)'^ are their Huang- Rhys factors-^^. From 
Eqs. Q and Q, the initial (fast) injection of an exci- 
ton, either coherently or incoherently, creates a sudden 
force on the modes, which we assume to be initially in a 
thermal, i.e. incoherent, state. This sudden force initi- 
ates transient oscillations of the modes at approximately 
their natural frequency ujk. To flrst order these transients 
can be treated as coherent oscillations - an assumption 
that will be corroborated later in our exact numerical 
treatment - and their back action on the excitons then 
acts essentially like a time-dependent fleld that drives co- 
herent, Rabi-like transitions between dissipative exciton 
states via the Hamiltonian terms 
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where {Qnm)it) oc ^ik^/%uJkCnCl^sm{uJkt) in the 
above-mentioned approximation of the initial, transient 
and coherent response of the modes to exciton injection. 
Heuristically, this non-equilibrium, laser-like driving es- 
sentially generates new electronic coherences in the sys- 
tem to replace those that are continuously damped out by 
the fluctuations of the smooth background environment. 
We note that in the language of chemical physics, this 
phenomenon arises through coherent non-adiabatic cou- 
pling (via discrete mode motion) which induces oscilla- 
tory crossing of potential surfaces. The actual (coupled) 
motion of the excitons and modes is more complex, and 
may also interact with non-Markovian dynamics of the 
background, but the physical picture presented here illus- 
trates a key point: electronic coherence may emerge from 
transiently exciting robust, weakly dephasing vibrational 
coherences which are then used to later back-transfer co- 
herence to the excitons, a novel type of coherence gener- 
ation and storage (se#2H^ for related observations that 
non-equilibrium systems may generate or maintain quan- 
tum entanglement). 

We should stress that this (re)generating of electronic 
coherence is a very different concept from protection 



of coherences by spatial correlations, non-Mar kovianity, 
or intrinsically weak environmental dephasing^^ ^^ ^^ ^^. 
Crucially, mode-driven coherences will be prominent 
whenever vibrational modes have frequencies compara- 
ble to exciton energy differences of strongly coupled 
chromophores and have dephasing times on picosecond 
timescales. This is of particular importance in PPCs, 
such as the FMO complex, since it allows regenera- 
tion of excitonic coherences whose energy splitting is 
commensurate with the frequencies of low energy envi- 
ronment modes. Many examples of such modes have 
been observed experimentally in PPCs, such as the 
FMO complex^^^, but their importance for interpret- 
ing experimental observations in multidimension al spec- 
troscopy has only just begun to be appreciate d^^*^^ * ^ -^. 
This mechanism, based on resonant driving of inter- 
exciton transitions, will not regenerate or sustain ground 
to excited state coherence, as the environment does not 
couple these states. This provides a natural understand- 
ing of the distinct coherence time scales in the problem, 
and we finally note that the tools and ideas we develop 
here may also be of relevance for photosynthetic charge 
transfer dynamics, where strong coupling between elec- 
tronic states due to coherent nuclear motion has been 
shown to enhance the efficiency of charge separation^^ . 

Fundamental mechanisms supporting long lived co- 
herence: Semiclassical results for the Fenna- Matthew s- 
Olson complex - To illustrate the physical mechanism 
of coherence regeneration through coherent (determinis- 
tic) motion of discrete modes in excitonic dynamics, we 
now present semiclassical simulations obtained by treat- 
ing the coordinates and momenta of the discrete vibra- 
tions as classical variables. This creates an effective time- 
dependent Hamiltonian for the excitons, and we use this 
in a secular Bloch-Redfield master equation with time- 
dependent dephasing rates for the reduced exciton state 
p(t), taking the smooth part of the experimentally- fitted 
spectral function of Adolphs and Renger (AR) for the 
background bath^^. In the following simulations, inter- 
actions with a mode of frequency 180cm~^ and Sk = 0.22 
were included, and the explicit form of the spectral func- 
tion is given in Appendix ] A 1[ The mode parameters and 
dephasing rate (lps~^) are taken from hole burning or 
estimated from fluorescence line narr owing ex periments, 
as is the background spectral functioiP^^^^l^. The back- 
ground bath and modes are initially in thermal equilib- 
rium at temperature T, the background reorganisation 
energy is A = 35 cm~^ and the seven-site FMO Hamil- 
tonian of C. tepidum was taken fromr^. We denote the 
matrix elements {en\p{t)\em) = Per.emi'^)- ^ ^^^^ descrip- 
tion of this semiclassical approach and the importance of 
time-dependent pure dephasing for super-Ohmic spectral 
densities is set out in Appendix [B] 

Figure [2] A shows results at T = 77K. Motivated by 
current experimental dataA-, which focuses on the two 
lowest exciton states, we present in Figure [5] (bottom) the 
coherences Re[peie2(^)] between exciton states (|ei), |e2)) 
with a symmetric superposition of these as the initial 
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FIG. 2: In all plots (A-D) results including and excluding the 
resonant 180cm~^ mode are shown as full and dashed lines, 
respectively. A&B. Semiclassical results at T = 77K. A. 
Inter-exciton coherence Re[peie2(^)] are shown for initial ex- 
citon state 2p(0) = |ei)(ei| + |e2)(e2| + |ei)(e2| + |e2)(ei|. B. 
Ground-excited coherences Ahs[peig{t)],Ahs[pe2g{t)] for ini- 
tial exciton state2p(0) = |e^)(e^| + \g){g\ + \ei){g\ + \g){ei\{i = 
1,2), respectively. Note that the optical high-frequency os- 
cillations of the peigit) coherences have been suppressed by 
taking the absolute value. Exciton population dynamics fol- 
lowing the injection of an excitation on site 1 of the FMO 
complex were also computed at this temperature, showing 
that excitons relax to the lowest energy states localised around 
sites 3 and 4 after ?^ 3 ps (not shown). This is inline with 
the experimental transport times of several picoseconds^^ . C 
& D. Semiclassical results at T = 277K. The same quantities 
and initial conditions as A &; B are plotted and used in C 
& D, respectively. The results at T = 277K in C &; D are 
plotted on the interval [0.2,0.8] ps to highlight the long lived 
coherences in the presence of the resonant 180cm~^ mode on 
timescale relevant for experiments at this temperature. 



state to resemble the laboratory condition after excita- 
tion with multiple laser pulses. Note that these results 
are obtained from a full simulation of the FMO complex, 
and so the influence of all sites on the lowest two exci- 
tons are included. In the presence of the 180cm~^ modes, 
which are nearly resonant with E1—E2 = 150cm~^, after 
a transient behaviour in the first 250fs the coherence ex- 
hibits prominent oscillations with an effective coherence 
time of ~ 950 fs. This is consistent with the coherence 
times of beating signals seen in FMO experiments^. In 
contrast, for the same background spectral density but 
neglecting the 180cm~^ mode, oscillations decay with 



a coherence time of just 280 fs after 250 fs. A multi- 
frequency beating pattern in the oscillations is appar- 
ent in the presence of resonant modes until about 850 
f. This is caused by the beating and dynamical cou- 
pling between the electronic coherent oscillations caused 
by the initial laser excitation and those induced later by 
the discrete mode motion. The vanishing of this pat- 
tern at the same rate as the exciton coherence in the 
background-only simulation provides direct evidence that 
initial coherences are not protected but replaced by mode 
interactions at later times. Both results also show a fast, 
time-dependent pure dephasing component at early times 
due to the background environment, described in detail 
in Appendix [B) 

Figure |2] B also shows the absolute value of the co- 
herence between |ei^2) and the optical ground state |^), 
starting from a symmetric superposition of |ei^2) and \g) 
which closely resembles the evolution probed through 
accumulated photon echo experiments^^ and recently 
developed single molecule femto-second pulse shaping 
techniques'^. An additional faster component to the co- 
herence decay is seen for both coherences at early times 
when the 180cm~^ mode is included, but the dynam- 
ics are qualitatively similar to the background-only case. 
The residual, slowly decaying component of pe^g is due 
to the relatively long (^ 2.4 ps) lifetime of the lowest en- 
ergy state at T = 77 K and the absence of pure dephasing 
in the long time limit in the AR spectral function in the 
Markov approximation - see Appendix |B 1| 

Figures |2] C and D show corresponding results at 
T = 277 K. Coherent oscillations lasting up to at 
least 600 fs have been observed in FMO^ at this near- 
ambient temperature. Here we find that the mode in- 
teractions lengthen the effective coherence time relative 
to the background-only simulations even more dramati- 
cally, with mode-induced coherence lasting up to 800 fs, 
compared with just 200 fs in the absence of the mode. 
The pei 29 coherences all decay with a similar, monoex- 
ponential time constants in the range 80 — 100 fs, which 
are dominated by the short (and almost equal) lifetimes 
and transient pure dephasing rates of the exciton state 
populations at T = 277K. 

Numerically- exact simulations- The essentially 
Markovian treatment presented so far has served to 
illustrate the principal physical mechanism responsible 
for long-lived electronic coherences represents an approx- 
imation whose validity must be assessed. Therefore we 
now present numerically exact results which include all 
possible effects of discrete mode motion, non-adiabatic 
coupling and fiuct nations, as well as the non-Mar kovian 
background. For clarity the new physics are presented 
for a dimer PP G, an important component of a range 
of natural PPCs^^^^. We again consider the same 
background spectral density J{uj) and 180cm~^ discrete 
mode of AR^. A 37cm-i ^^^^ ^^-^^^ g^ ^ ^ l^^ ^j^-^j^ 

has recently been used to_ describe features of 2D FMO 
spectra is also included'^, with parameters taken from 
fiuorescence line narrowing experiments^^. Electronic 



parameters used in the simulations are J12 = 53.5cm~^, 
61 — 62 = 130cm~^. This gives two exciton eigen- 
states with an energy difference of 170cm~^, which 
is based on sites 3 and 4 of the FMO HamiltoniarP^, 
and which are fairly typical values for exciton pairs 
in PPCP^El. The evolution of the global system- 
environment density matrix is then computed using 
a new finite-temperature extension of the numerically 
exact TEDOPA methocPK^- see Appendix |A2J In all 
simulations the initial state is a product of an exciton 
state and a thermal state of the environment (which 
includes the discrete modes). To show that the main 
effects arise from the presence of the resonant 180cm~^ 
mode, simulations are carried out both with and without 
this mode, whilst the non-resonant 37cm~^ mode is 
always retained. 

Inter-exciton and ground- exciton coherences- Figure [3] 
A shows the evolution of the electronic coherence peie^ if) 
starting from an initially prepared (pure) symmetric su- 
perposition of the exciton states |ei) and |e2) at T = 77 
K. With coupling to the 180cm~^ mode, the multi- 
frequency beating and revival dynamics in the oscillating 
coherence are again seen, indicating mode-driven coher- 
ence. This leads to coherence oscillations with a fast ini- 
tial dephasing time (200 fs) and a residual, longer lasting 
component arising from mode driving. Similar revival 
patterns have been observed in time resolved spectra of 
FMC0^. At T = 277 K, shown in Fig. [3]C, the resonant 
mode also greatly enhances coherent oscillations relative 
to the background-only coherences over the first 800 fs, 
which is, again, consistent with experimental results^. 
The non-resonant mode at 37cm~^ plays no significant 
role in these inter-exciton dynamics, as expected. 

Figure [3p shows the absolute value of the ground- 
excited state coherences pcngf) at T = 77 K for ini- 
tial symmetric superpositions of states \en) and \g). The 
thermal and quantum fiuctuations of both the 180cm~^ 
and 37cm~^ modes induce amplitude modulations of 
these coherences at the mode frequencies'^, leading to 
an effectively faster and oscillatory initial decay and then 
slow oscillations on top of the decaying amplitudes. The 
short effective time constants for all curves lie in the 
range 150—200 fs, which are consistent with experimental 
results , and were faster than those found in the semi- 
classical approach. The coherence pe2gf) decays through 
population relaxation from |e2) to |ei), however, on top 
of the periodic modulation caused by modes, additional 
weak recurrence features result from oscillatory popula- 
tion transfer induced by the 180cm~^ mode - see Fig. 
p] The coherence peig{t) shows a fast, incomplete and 
non- markovian decay (~ 150 fs) initially, and then an 
extremely slow decay due to the long-time absence of 
pure dephasing in the AR (super-Ohmic) spectral den- 
sity and the long (^2.4 ps) lifetime of the lowest energy 
excited state |ei) at T = 77 K. Long \ei){g\ coherence 
times (10 — 100 ps) have been seen in previous FMO 
photon echo measurements at low temperatures (T < 50 
K)^^ At T = 277 K (Fig. [3]D), both ground-excited 
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FIG. 3: In all plots (A-D) results including and exclud- 
ing the resonant 180cm~^ mode are shown as full and dashed 
lines, respectively. A&B. TEDOPA results at T = 77K. A. 
Inter-exciton coherence Re[peie2(^)] are shown for initial ex- 
citon state 2p(0) = |ei)(ei| + |e2)(e2| + |ei)(e2| + |e2)(ei|. B. 
Ground-excited coherences Ahs[peig{t)],Ahs[pe2g{t)] for ini- 
tial exciton state2p(0) = |e^)(e^| + \g){g\ + \ei){g\ + \g){ei\{i = 
1,2), respectively. Note that the optical high-frequency os- 
cillations of the Peigit) coherences have been suppressed by 
taking the absolute value. Weak revivals of Abs[pe2p(^)] (o^ 
top of faster amplitude modulations) in the interval [0.3, 0.6] 
ps arise from coherent population transfer {against the en- 
ergy gradient) and match similar features seen in Fig. HI A. 
C &; D. TEDOPA results at T = 277K. The same quantities 
and initial conditions as A &; B are plotted and used in C 
& D, respectively. The results at T = 277K in C &; D are 
plotted on the interval [0.2,0.8] ps to highlight the long lived 
coherences in the presence of the resonant 180cm~^ mode on 
timescale relevant for experiments at this temperature. 



coherences decay rapidly due to the fast population re- 
laxation to thermal equilibrium and enhanced pure de- 
phasing. Again, the additional oscillatory decay in the 
presence of the resonant modes in the TEDOPA simula- 
tions leads to significantly faster reduction of coherence 
(with time constant of 50 fs) compared to the semiclassi- 
cal approach which neglects the quantum fluctuations of 
the environment. 

Population oscillations and environmental variables- 
For T = 77K, Fig. |4] A shows the population remaining 
in 1 62) for the initial state p(0) = |e2)(e2|, which does not 
contain any initial inter-exciton coherence and may be 
created by incoherent excitation alone. The presence of 



the 180cm~^ mode induces population oscillations in the 
population decay, which explain the similar oscillations 
seen in peig{t) in Fig. [sIB, and is also direct evidence 
of weak coherence transfer (not shown). Population os- 
cillations have also recently been observed in FMCP^, 
and imply a partially reversible energy exchange with 
the "environment" that transiently violates detailed bal- 
ance. Figure [4] B shows the spontaneous generation of 
oscillatory coherences over these population dynamics. 
Non-Markovian effects (finite reorganisation time) and 
the 37cm~^ mode are also seen to generate short-lived 
spontaneous coherences, however, these are transient and 
only the 180cm~^ mode drives strong long-lasting coher- 
ence oscillations. No spontaneous coherence arises in the 
standard Bloch-Redfield description, showing the essen- 
tial need for going beyond such theories - See Appendix 
B 1[ In Figure |4] C, we make use of the complete envi- 



ronment information provided by TEDOPA to show the 
expectation value of the collective environment coordi- 
nate Xi^^^EU. Figure [4] C shows that the exciton-mode 
interaction leads to long-lasting coherent environmen- 
tal oscillations which drive exciton transitions, confirm- 
ing the semiclassical picture and coherence-generating 
mechanism previously outlined. In the absence of res- 
onant modes, these strong oscillatory features vanish in 
the collective mode and exciton populations. Figures [4] 
D-F show that, remarkably, the same non-equilibirum 
and spontaneously coherent dynamics are also present 
(though less strong) at T = 277K. 

We finally observe that the spontaneous coherence in 
Fig. ffl B tends to a (small), non-oscillating, non-zero 
value at long times at T = 77 K, indicating that the 
exciton-mode steady states generated by the dynamics 
are quantum superpositions of excitonic and vibrational 
degrees of freedom. This residual coherence vanishes at 
T = 277 K. More striking examples of how mode in- 
teractions may alter relaxed electronic states is given in 
Appendix [C| which sets out how discrete modes may pro- 
vide insight into anomalous oscillations recently observed 
in conjugated polymerd^^^. 

Single system vs ensemble- In the linear absorption 
spectra of FMO ensembles, the lowest absorption band 
- thought to consist of just a single excitonic transition, 
and possibly some vibrational sidebands^^ ^^ ^^ - has an 
inhomogenous line width of 80cm~^ FWHM^^. This im- 
plies an inhomogeneous dephasing time of the \ei){g\ co- 
herences of 130 fs, which is close to the experimental 
measurements at 77 Kp'. In ensemble 2D spectroscopy 
experiments, this might completely mask the long (sin- 
gle complex) coherence time of |ei) (^| at T = 77 K in Fig. 
|3]B, and - moreover - if the inhomogeneous broadening 
of site energies on each pigment is independent, then sig- 
nals from inter-exciton coherences such as |ei)(e2| might 
also dephase on similar timescales due to the resulting 
width of inter-exciton energy differences^. 

The actual observation of long-lasting inter-exciton 
coherences in ensemble experiments thus requires fur- 
ther explanation, and could be supported by correlated 
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FIG. 4: A. Population dynamics, showing pe2e2{t) for elec- 
tronic and mode parameters in the text. Initial state was 
p(0) = |e2)(e2| and environment was at T = 77K. Weak, oscil- 
latory revival of population in the |e2) state between 0.3 — 0.6 
ps matches the revival dynamics of pe2g in Fig- [3]B, indicat- 
ing coherent population (and coherence) transfer. B. spon- 
taneous electronic coherences Re[peie2{t)] for same parame- 
ters and initial conditions. C. Average displacement (Ai)(t) 
for same parameters an initial conditions. D-F. The same 
as A-C, respectively, but with the environment initially at 
T = 277K. 



and incoherent initial conditions. In focusing on purely 
electronic observables - which are the relevant degrees 
of freedom probed in actual EET experiments - we have 
traced over the quantum fine structure of the hybrised 
exciton-mode states, and have shown that an effective 
semiclassical motion of quasiresonant modes already pro- 
vides a qualitative and intuitive solution to the problem 
of how both long and short lasting electronic coherences 
may naturally appear in the PPCs. This heuristic model 
is subsequently corroborated by means of an exact nu- 
merical treatment that includes naturally the quantum 
fluctuations in the vibrational environment. This insight 
illustrates a situation outside the standard open quan- 
tum system paradigm, where parts of the environment 
(the sharp vibrational modes) are driven strongly out of 
equilibrium and counteract effects induced by the ther- 
mal background. The matching and coupling of under- 
damped vibrational frequencies to excitonic transitions 
presents a clear new design principle for building coher- 
ent exciton transport structures, and thus raises new pos- 
sibilities within the context of noise- assisted transport, 
such as engineering discrete bath structures for (coher- 
ent) energy transport -recently exemplified by the con- 
cept of the phonon antenna and recent studies of en- 
hanced non-equilibri um en ergy harvesting efficiency away 
from detailed balance^^J^. From the point of view of fun- 
damental physics, the experimental verification of the dy- 
namics discussed in this work would provide additional 
ground for the consolidation of quantum biology as a 
new, truly multidisciplinary research field. 
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Appendix A: Spectral functions, excitonic 

parameters and the finite-temperature TEDOPA 

method 



1. Spectral function 

The spectral function used to characterise the environ- 
ment acting on each site in the TEDOPA simulations is 



given by 

2 

J{uj) = Jo{uj)^'^SiUj^S{uj-uJi), (Al) 



Jo{uj) 



X [lOOOcj^e" 



4.3uj^e~ 



9!(1000cj5+4.3cjg) 



(A2) 



The bath parameters used throughout are Si = 
0.12, ^2 = 0.22, cja = 0.57cm-i ^^^ ^^ ^ 1.9cm-^ 
T hese parameters were extracted from experiments 
^j j34 | 43 [ Vibrational modes with similar frequencies have 
a lso be en experimentally observed in the FMO complex 
^j ji8|2Q [ rpj^^ reorganisation energy used in our simula- 
tions is A = 35cm~^, taken again from^^. The same 
spectral function is used for the semiclassical simulations, 
except the 37cm~^ mode is excluded. 



2. Simulation parameters and method 

The parameters of the dimer system used through- 
out are 61 — 62 = 130cm~^, and J = 53.5cm~^. Di- 
agonalising the dimer Hamiltonian yields two singly- 
excited exciton states |ei,2) with an energy splitting 
£^12 = E1—E2 = 170cm~^. Each constituent of the dimer 
is linearly coupled to an environment of harmonic oscil- 
lators as in the well-known spin-boson model. This prob- 
lem is then transformed to a physically equivalent setting 
where each constituent of the dimer is coupled to a lin- 
ear chain of harmonic oscillators with nearest neighbor 
interactiorP^Hill^ The initial state of the global system 
is always taken as ptotai = p(0) ^ ^-^B/kkT -^^ ^ prod- 
uct of an arbitrary exciton state and the thermal state 
of the uncoupled environment at temperature T. The 
initial thermal environment is prepared using imaginary 
time evolution and the coupling between the system and 
thermal environment is turned on at time t when the ini- 
tial thermal state of the environment has been prepared. 
Hence the problem is now the simulation of the time evo- 
lution of this 1-D quantum system. To this end the global 
state of the system and environment is then evolved using 
a finite-temperature t-DMRG algorithm which retains in- 
formation of the state of system and environment by com- 
bining matrix product states (MPS) with a sophisticated 
way of implementing local updates to the system yield- 
ing an efficient simulation of one-dimensional quantum 
systems (see Zwolak et aP^ for the description of imag- 
inary time evolutions as well as real time evolutions in 
this framework). The time dependence of arbitrary ob- 
servables O of the excitons or bath are then obtained 
from {0){t) = Tr[p(t)0], where p{t) is the total state of 
the system and the bath. For the present simulations 
the chains representing the environment are 49 sites long 
(preventing recurrence effects on the relevant time scale) 
and are composed of harmonic oscillators truncated to 
three levels. Using a bond dimension of 320 for the MPS 
covers all essential physical effects and restricts the re- 
quired numerical resources to a reasonable limit. The 



simulations have been carried out using a second order 
Suzuki- Trotter expansion with time steps of 1/188 and 
1/136 of the total time for the imaginary and real time 
evolution, respectively. These parameters have been ad- 
justed by observation of the norm during the simulation, 
as to keep the incurred errors due to the truncation of 
the MPS and the Suzuki- Trotter expansion to a tolera- 
ble value. A detailed account - with several applications 
- of this new finite-temperature version of the TEDOPA 
approach to open systems will be presented elsewhere. 



3. Exciton basis Hamiltonian and equations of 
motion 

Rewriting the Hamiltonian of Eq. (IT]) in the exciton 
basis \i) = Xl^C'^len), we obtain, 

H = ^En\en){en\ +y^Qnm|en)(em| 
n n,m 

+ ^ujkalf^aik. (A3) 

ik 

The new couplings to the bath are given by Qnm = 
^i,k 9ikCl,Cl^{aik + a]^). The coefficients C'^ can be 
taken as real, so Qnm = Qmn- All bath modes are 
taken to be independent, i.e. [ai/c,aU = 6ijSkj^ where 
the square bracket represents a commutator. In this ba- 
sis, we see that the environment now couples to both 
the energy level populations (diagonal coupling) and can 
induce transitions between exciton states (off-diagonal 
coupling). Although transverse terms are important for 
the semiclassical driving of new exciton coherences by 
the modes, the diagonal terms are important in the ini- 
tial excitation of the mode oscillations, especially when 
starting from an exciton population state. 



Appendix B: Semiclassical equations of motion 

We have developed a theory in which the key process 
mediated by the discrete modes can be interpreted as a 
quasi-coherent driving of the exciton system by the co- 
herent displacement of the modes. This suggests that 
a significant part of the interactions between the dis- 
crete modes and the excitons could be accounted for in 
a semi-classical approximation. To derive the semiclas- 
sical equations of motion which were used to investigate 
exciton dynamics in the FMO complex, we treat the dis- 
crete modes as part of the system and account for the 
effects of the background environment using a simple 
Bloch-Redfield approach. The effective system Hamil- 
tonian is the same as Eq. (A3), but the sum over bath 



modes is restricted to just the two discrete modes that 
appear in the spectral function of Eq. ( |A1[ ) . We then use 
the Heisenberg picture to compute exact equations of mo- 
tion for the mode displacements Xik{t) = {ciik{t)-\-(\k{'^)) 
and their momenta Pik{t) = {(iik{t) — (^Iki'^))^ where 



A{t) = e^^t Ae~^^^ denotes operators in the Heisenberg 
picture. The expectation value refers to the operation 
{A) = TT[ptotai{t)^]- In this notation the index i refers 
to modes coupled to site i and the k index labels the 
different discrete modes at each site i.e. the modes of 



frequency ujk in Eq. (Al). As an example, the equation 



of motion for the momentum operator is 



dt 



{^^k(t) - a\^(t)) = i[H, {a,k{t) - aim (Bl) 

= -iuji{aik{t) ^al{t)) 

+ 2i^/S~kOOk^\Ci,\^\en){enm 



+ 2zVSkUJk}^C:,Cl\em){enm. 

Taking the expectation value of this equation of motion 
using the result (|e,n)(en|(^)) = T^Aptotai{t)\em){en\] = 
Pe^ 6^(^)5 where pe^e^(t) are the matrix elements of the 
reduced density matrix of the excitons, we arrive at the 
equation of motion for the expectation values 



dPik{t) 
dt 



-iuJiXik{t) 



Pe 



At) 



2iy^k00k Yl ^nQ 



,Pe 



n^m 



and analogously for the position operator 

dXik{t) 



dt 



-iuJkPik(t)' 



M (B2) 



(B3) 



The semiclassical approxi mati on amounts to replacing 
the mode operators in Eq. ( A3 ) by these time-dependent 



c- number expectation values, creating a time-dependent 
Hamiltonian Hxi{t) for the excitons. We evolve the re- 
duced density matrix of the exciton system according to 



dpsjt) 
dt 



-i[HxM.pm^nt)[ps{t)]. (B4) 



where we have added a dissipative term Tl{t) that ac- 
counts for a background environment described by the 
smooth background part of the spectral function of Ap- 
pendix A 1 The Redfield relaxation tensor 7l{t) in Eq. 



(B4) is computed in the absence of the discrete modes 



and is restricted to secular (energy-conserving) terms 
only. The full equation and expressions for the relaxation 
and d ephas ing rates contained in 7l{t) are given explic- 
itly i] j^5 | 46 [ However, the time argument is included in 
Tl{t) to denote that some of the rates - corresponding 
to pure dephasing - are non-standard, and are actually 
time-dependent. This will described in Appendix |B 1[ 

In all simulations the initial mode displacements and 
momenta are taken drawn from a thermal distribution 
at temperature T. A phenomenological damping rate 
7 for the modes, which is assumed to be the same for 



both modes, is set to g ive a typical da mpin g time for 
vibrational motion of 1 p ji'^ l 34[ xhus Eq. (B2 ) is replaced 

by 



dt 



-2'yPik(t) -iuJiXik{t) 

2iV^C^fc^|C;|Ve„e„(t) 
n 

2i^/S~kC0kY.CnClnPe^eAt)■ (B5) 



For completeness, the equation of motion for the density 
matrix elements pcr^emit) = {^n\ps{t)\em) are 



dper^emjt) 
dt 



-iEnmPer,emit) " ^ ^^ Qnj ft)Pe,e^ ft) 
J 

-\-i ^ QjmPer^ej (t) + ^ ^nmjkPejCi (^), (B6) 



where Qnm{t) = Y.i^k^/^^kClCl^Xik{t). From Eq. 
( |B6[ ) we see that the mode displacements lead to time- 
dependent couplings of populations and coherences which 
generate the 'driving' of inter-exciton transitions which 
are shown in the main text to enhance coherence lifetimes 
and generate new coherences. Oscillatory population dy- 
namics also appear through the coupling of populations 
to coherences mediated by the discrete modes. One can 
immediately see from the structure of these equations 
that discrete modes cannot drive \en){g\ coherences, and 
thus do not regenerate or support ground-excited coher- 
ence. However, these equations do predict pe^g -^ Pcmg 
coherence transfers due to coherent motion of the discrete 
modes which lead to rapid suppression of coherence ex- 
cept for the lowest exciton. For each set condition of 
initial mode conditions drawn from the thermal distribu- 
tion, the equations of motion for Xi{t)^Pi{t) eqs. |B3|B5 
and eq. |B6| are integrated together, and these trajectories 
are averaged over to obtain the results of Fig. [2] 

We finally note that in this simple description, quan- 
tum fiuctuations of the discrete modes are neglected, as 
are the renormalisation effects of the background envi- 
ronment. This last approximation is appropriate for the 
relatively weak coupling parameters used for the back- 
ground spectral density. All of these effects are captured 
in the TEDOPA simulations. The main qualitative dif- 
ferences were observed in the \en){g\ coherences, where 
amplitude modulations due to quantum fiuctuations of 
the discrete modes could not be reproduced by the semi- 
classical approach, leading to an overestimate of coher- 
ence times. 



1. Non-Markovian Pure dephasing rates for 
super-Ohmic spectral densities 

In standard, secular Bloch-Redfield theory, the Marko- 
vian pure dephasing rate Tpd (which causes decoherence 




T=77K ■ 
T=277K ■ 



Time / ps 



FIG. 5: pi2(t) =p(0)e-^orpd(t')dt'_ initial condition is p(0) = 
1 and the transient, non-Mar kovian pure dephasing rate is 
computed using Eq. (B8) for the spectral density given in 



[AT] Results are shown at T = 77K (black line) and T = 277K 
(blue line) . The vanishing of the dephasing rate at long times 
leads to the survival of a finite value of coherence at long 
times. 




FIG. 6: Abs[pei^(t)], for the standard, secular Bloch-Redfield 
theory (Black line) and for the generalised theory which in- 
cludes a time-dependent pure dephasing rate (Blue line). Ini- 
tial condition is 2y9(0) = \ei){ei\ + \g){g\ + \ei){g\ + \g){ei\ and 
the bath is at T = 77K. Calculations are for the fuh FMO 
system described in the main text. 



of energy eigenstate superpositions in the absence of en- 
ergy relaxation) at temperature T for a spectral density 
J{uj) is (neglecting matrix elements) proportional to 



- pd 



oc lim 



dt 



oc lim 



duoJiuo) coth 



puj 



puj 



COS {{uo — e)t) , 



/>oo 

im TT / dujJ{(jj) coth ( ^ ) 5{uj — e), 



oc lim 



OC lim 



TT J(e) coth 



2jj(e) 



/?e 



(B7) 



where j3 = (ki^T) ^. The pure dephasing rate is thus con- 
trolled by the behaviour of spectral function approach- 



ing zero frequency^ -^^^', as one would expect for a pro- 
cess that involves no energy dissipation. Typical spectral 
densities have power-law dependence at low frequency 
J{uj) oc o;^, with s some exponent. For super Ohmic 
damping (s > 1) the limit in Eq. (B7) is zero, and 



thus pure dephasing is apparently absent for this class of 
spectral functions^^. This can considerably weaken the 
temperature dependence of dephasing, which can now 
only occur through energy relaxation in Bloc h-Redfield 
theor}/ — . The AR spectral function of A 1 is in this 
class. The commonly encountered Ohmic spectral den- 
sity (5 = 1) leads to a finite limit and a pure dephasing 
rate proportional to temperature^^. The relevance of this 
observation h as rece ntly been discussed in the context of 
PPG dynamics^^ ^^, and the absence of pure dephasing in 
super-Ohmic spectral densities has been put forward as 
a means of prolonging excitonic coherence whilst main- 
taining efficient energy transport^^. 

However, the pure dephasing rate goes to zero in the 
full Markov approximation due to the fa ct t hat the up- 
per limit in the integral over time in Eq. (B7) is taken to 



infinity, leading to the delta function that picks out the 
zero frequency component. However, a more accurate 
treatment would take the upper limit to be t (finite) so 
that the (now time-dependent) pure dephasing rate also 
remains finite (it samples a range of finite frequencies 
proportional to 1/t) at short times, and, at least at those 
early times, may contribute significantly to the dephas- 
ing of coherences. To take account of this, we therefore 
introduce a time-dependent pure dephasing rate Tpd{t) 
in place of the vanishing Markovian rates in our Redfield 
tensor proportional to 



^pd{t) (X dt dujJ{uj) coth [^) cos{ujt). (B8) 

The key point is that while Tp^it) ^ as t ^ 00, the 
evolution of a coherence pi2{t) obeying the equation of 
motion pi2{t) = -Tpd{t)pi2{t) iF^ 



Pi2{t) = Pi2{0)e-fo^^^it')dt'^ 



(B9) 



and the exponential suppression resulting from the inte- 
gral over the transient time-dependent rate can be non- 
negligible. Figure [5] illustrates this, showing this sup- 
pression factor for the AR background spectral density 
at T = 77 K and T = 277K, taking Tpd{t) to be equal to 
the RHS of Eq. (B8|. Theses curves shows a fast initial 
suppression which then stops and plateaus as the pure 
dephasing rate vanishes. The 'residual' coherence is sup- 
pressed at higher temperature, in line with expectations 
based on exactly solvable, non-Mar kovian pure dephas- 
ing models which also predict the residual coherence for 
super-Ohmic spectral densities^. Despite the long-time 
absence of pure-dephasing. Fig. |5] clearly demonstrates 
the importance of accounting for the early time suppres- 
sion, and naive applications of standard Bloch-Redfield 
theory for super-Ohmic spectral densities are likely to 
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FIG. 7: Ahs[peig{t)], for the standard, secular Bloch-Redfield 
theory (Black line) and for the generalised theory which in- 
cludes a time-dependent pure dephasing rate (Blue line). Ini- 
tial condition is 2p{0) = \ei){ei\ + \g){g\ + \ei){g\ + \g){ei\ and 
the bath is at T = 77K. Calculations are for the full FMO 
system described in the main text. 



over-estimate effective coherence times without this cor- 
rection. 

We illustrate this further with Fig. [7| which shows 
Abs[pe-^^g{t) for the full FMO complex computed with 
and without the super- Ohmic time-dependent pure de- 
phasing. The Bloch-Redfield lifetime of the (lowest ex- 
citon energy) |ei) state at T = 77i^ is 2.4 ps, imply- 
ing a very long dephasing time of 4.8 ps in the absence 
of pure dephasing. This is completely inconsistent with 
the TEDOPA results in Fig. |3J while the corrected non- 
Markovian dephasing result actually gives much better 
agreement w.r.t overall dephasing rates and the long-time 
value of the coherence. One sees that the dephasing is 
dominated by the early-time non-Markovian correction, 
with the much weaker dephasing arising from the relax- 
ation of the |ei) state leading to a slow decay of the 
residual coherence plateau. 

We note, however, that the inclusion of a time- 
dependent, non-mar kovian pure dephasing rate in our 
simple time-local master equation is an approximation, 
and for some spectral densities the time-local master 
equation may be invalided. Fortunately, this important, 
transient pure dephasing effect is handled exactly by 
TEDOPA, and the results of this section are provided 
to show the necessity of using advanced techniques for 
the open-quantum dynamics of PPCs, as well as the im- 
provements to numerically cheaper techniques that in- 
sight from TEDOPA can generate. 



Appendix C: Anomalous anisotropy decay in 

organic polymers and spontaneous generation of 

coherences 



In a recent paper Collini and Scholes demonstrated 
room-temperature inter-exciton coherence in conjugated 
polymers in chloroform solutiorP^^. They present some 




Time/ ps 



FIG. 8: Exciton 2 population for initial condition ps{0) = 
|e2)(e2| using the semiclassical equations at T = 77K, tyib = 1 
ps and S = 0.00,0.22,0.50,1.00 (black, red, blue, magenta, 
respectively). Horizontal axis denotes time in ps. 



anomalous anisotropy decay data in which long-lasting 
(> 1 ps) coherent oscillations and bi-exponential popu- 
lation transfer in the anisotropy decay were observecP^. 
These oscillations had frequencies consistent with vibra- 
tional frequencies in the system, but were also rather 
similar to exciton energy differences as well. The data 
showed conflicting features (lack of coherence time de- 
pendence, long-lasting oscillations) which meant that it 
was not possible to assign these features unambiguously 
to either electronic or vibrational coherences. The model 
and mechanisms presented here may provide some in- 
sight into the these data, as they show how the gen- 
eration of spontaneous coherences by the mode-exciton 
interaction may induce population oscillations on long 
time scales even when the exciton system is populated 
without exciting vibrational coherence. Figure [8] shows 
some illustrative population dynamics obtained using the 
semiclassical approach of section [S] for a dimer with 
61—62 = lOOcm"^, Ji2 = 86cm~ , each site of which 
is coupled to a discrete mode with uji = 180cm~^ and a 
variable Huang-Rhys factor S. The background bath is 
at zero temperature. Two key features which emerge in 
the presence of the mode interaction is that the popu- 
lation decay develops a number of new timescales (rel- 
ative to the case Si = 0.0). The mode leads to a 
multi-exponential population decay (initially faster when 
modes are present) on which complex coherent oscilla- 
tions resulting from the mode interactions are superim- 
posed. These coherent oscillations persist for at least the 
dephasing time of the modes (1 ps in these simulations), 
and are stronger and persist longer at stronger coupling. 
Figure [9] shows the spontaneous generation of the elec- 
tronic coherence by the discrete modes for the same 
parameters and initial conditions as in Fig. [8J This 
provides direct evidence that the population oscillations 
shown in Fig. p] are the result of inter-exciton (elec- 
tronic) coherences induced by the driving of the dis- 
crete modes. The magnitude, oscillation lifetime and fi- 
nal value of these spontaneous coherences increase with 
coupling strength, and show interesting dynamical os- 




Time / ps 

FIG. 9: Spontaneously generated coherences Re[peie2(^)] for 
initial condition ps{0) — |ei)(ei| using the semiclassical equa- 
tions at T = 77K, Uib = 1 ps and 5' = 0.00,0.22,0.50,1.00 
(black, red, blue, magenta curves respectively). Horizontal 
axis denotes time in ps. 
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dilations at very strong coupling, where weak signs of 
non-rotating wave harmonics can be seen in the driven 
coherences. The final values show that the discrete modes 
change the nature of the equilibrium exciton states and 
both exciton states are significantly hybridised by the 
modes at strong coupling. Note that these calculation 
are intended to illustrate the physics which might explain 
the data of Refs.^^ ^^, but the parameters we have chosen 
are fairly arbitrary and a much more detailed analysis of 
the physical and experimental set up is required to unam- 
biguously assign a mechanism to the observed behaviour. 
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